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ABSTRACT 



We investigate the extinction together with the radial velocity dispersion and 
distribution of red clump stars in the anti-center direction using spectra obtained 
with Hectospec on the MMT. We find that extinction peaks at Galactocentric 
radii of about 9.5 and 12.5 kpc, right in front of the locations of the Perseus and 
Outer arms and in line with the relative position of dust and stars in external 
spiral galaxies. The radial velocity dispersion peaks around 10 kpc, which co- 
incides with the location of the Perseus arm, yields an estimated arm-interarm 
density contrast of 1.3-1.5 and is in agreement with previous studies. Finally, 
we discover that the radial velocity distribution bifurcates around 10-11 kpc into 
two peaks at +27kms~ 1 and -4kms _1 . This seems to be naturally explained by 
the presence of the outer Lindblad resonance of the Galactic bar, but further ob- 
servations will be needed to understand if the corotation resonance of the spirals 
arms also plays a role. 

Subject headings: Galaxies: individual (Milky Way) — Galaxy: disk — Galaxy: 
kinematics and dynamics — Galaxy: structure — Stars: kinematics and dynam- 
ics 
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Introduction 



Non-axisymmetric rotating substructures, such as a bar and spiral arms in the Milky 
Way, induce a corotation resonance (CR) as well as i nner (ILR) and outer Lindblad reso- 



nanc e s (OLR) around cert ain Galactocentric radii (e.g.. lWeinberg 



19941 ; iFux |200(A iQuillen et al 



201 ll : lAntoja et al. 1 1201 if ). These resonances seem to be able to explain the moving groups 
in the velocity distribution in the Sola r neighborhood discovered from Hipparcos data and 



the Geneva-Copenhagen S urvey (e.g.. iDehnen I I200CH : iDe Simone. Wu. fc Tremaine 1 12004 



Guillen fc Minchev 1 120051 ; ISellwood 1 120091 ) . However, resonances from multiple patterns 
may exist and overlap, and moreover the Sun may not be located right on any of the res- 
onance radii. Henceforth, it is difficult to distinguish the effects of resonances in the Solar 
neighborhood. Observations of the velocity distribution outside the Solar radius are needed 
to constrain the resonances. 

We report in this letter a resonance feature at a Galactocentric radius of 10-1 lkpc in 
the radial stellar velocity distribution derived from MMT/Hectospec observations of 698 red 
clump stars in the anti-center direction. In Section [2] we briefly describe the observations, 
data reduction and analysis. In Section [3] we identify the location of the Perseus arm and 
discuss possible reasons for the bifurcation seen in the radial velocity distribution. Finally, 
we summarize our findings in Section HJ 



2. Data reduction and analysis 

We used MMT/Hectospec to obtain spectra with resolution R=2500 of stars in 4 one- 
degree-diameter fields between I = 178°-182° at b = 0°. Red clump stars are selected as 
targets from the 2MASS catalog by a single color cut J — K > 0.65 mag. Close-by dwarf 
stars can overlap in intrinsic color with the red clump stars, but as the latter get reddened by 
extinction in the mid-plane the dwarf stars can be easily removed. The single-color selection 
of red clump stars is thus expected to be very efficient with only little contamination by 
other giant stars and remaining dwarf stars. 

In total we obtained 993 spectra with mean S/N ~ 17. The reduced spectra, of which 
three examples are shown in the top panel of Fig. [TJ show that later- type K/M giant stars 
and dwarf stars contain stronger MgH and TiO absorption features within 5000A -5300A 
and 6000A-6500A than red clump stars. After computing the equivalent width (EW) of 
the absorption features in each spectrum, we arrive at the distribution in the bottom panel 
of Fig. [TJ with a Gaussian peak of most, if not all, red clump stars plus a tail of higher 
EW caused by the MgH and TiO bands in the contaminating giant and dwarf stars. After 
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Fig. 1. — Top panel: Reduced MMT/Hectospec spectra of a red clump star, dwarf star, and 
giant star. The equivalent width between the pairs of vertical dashed lines at 5000A-5300A 
and 6000A-6500A is used to remove the dwarf and giant star contaminants based on their 
prominent MgH and TiO bands. Bottom panel: The solid curve shows the distribution 
of the EW for all spectra. The dashed curve is the Gaussian fit to the peak with the 
dot-dashed vertical line indicating one standard deviation above the mean, above which 
stars are excluded from the red clump sample to safely remove all the dwarf and giant star 
contaminants in the tail. 



conservatively selecting the stars with E W not larger than one standard deviation above the 
mean, we are left with 747 red clump stars. This implies a completeness of at least 85% and 
hence confirming the high effectiveness of the color selection. To enable the derivation of 
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a smooth radial distribution in extinction and velocity below, we exclude another 46 stars 
with excessive local extinction as well as 3 stars that are 3-cr outliers in velocity, leaving a 
final sample of 698 red clump s tars. Three inde pendent methods, H a line fitting, synthetic 
template matching and ULySS (jWu et al. Il201ll ). are used to extract a mean radial velocity 
value and corresponding error from each spectrum. 



Next, we use the Hipparcos (jPerryman et al.lll997l ) and 2MASS ( iSkrutskie et al. 1120031 ) 
catalog to infer, after taking into account photometric and parallax errors, the intrinsic 
J — K color index and absolute i^-band magnitude of the red clump stars: (J — K)q = 
0.654 ± 0.058 mag an d Mk = —1.547 ± 0.311 mag . The latter is in agreement with both 
Groenewegen I (120081 ) and lLopez-Corredoira et al. ( 120021 ). and since metallicity is at most 
weakly correlated with Mj ( Zhao. Qiu &: Maol2001 ). and likely less so for M K , we can obtain 
accurate distance estimates to the red clump stars in the following way. 

For each red clump star the probability density function (PDF) of extinction given its 
observed color index, p(Ak\J — K), is computed, assuming it is a Gaussian distribution 
with dispersion given by the previously derived intrinsic J — K dispersion plus the J — 
K photometric error added in quadrature. Similarly, the likelihood of apparent .fT-band 
magnitude K given the distance modulus DM and Ak, p(K\DM, Ak), is assumed to follow a 
Gaussian distribution with dispersion given by the previously derived intrinsic Mk dispersion 
together with the photometric uncertainty in K. The PDF of the distance modulus given 
the observed J — K color and K then follows from Bayesian's theorem as 



p(DM\K,J - K) = Jp(K\DM,Ak)p(A k \J - K) p{A K ) p(DM) dA K , 



(1) 



with adopted uniform priors P(Ak) and P(DM). Adopting as Solar radius R Q 
distance modulus DM is converted to Galactocentric radius Rgc- 



8kpc, the 



Fig. [2] shows the distributions of K-band extinction Ak, radial velocity dispersions a^y 
and radial velocity (RV) versus Rgc distribution. Given a distance, the distribution in Ak 
follows by summing up p(Ak\ J — K)p(DM\K, J — K) for all stars. The line in the top 
panel in Fig. |5] shows the Ak with maximum probability. And the shadow shows the 68% 
confidence region. There are two significant jumps and hence dense dust regions around 
9.5 kpc and 12.5 kpc is found in the distribution in Ak- 

The radial profile of a^v in the middle panel is obtained from the weighted variance of 
RV, which equals to a RV plus the squared measurement error of RV, e 2 RV . The weight is given 
by p(DM\K, J — K) for each star with normalization. During the calculation we assume 
that the probability of the measured RV is a Gaussian with sigma e RV . The uncertainty of 
a^v is given from a bootstrap resampling process. It reveals that there is a clear bump at 
10 kpc. 
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Fig. 2. — Structural and velocity distribution of red clump stars as a function of Galacto- 
centric radius Rcc with the location of the Perseus and Outer arms indicated. Top panel: 
Most likely if-band extinction A K and its 68% confidence region. Middle panel: Same 
for the radial velocity dispersion o~rv- Bottom panel: The probability of radial velocity 
versus distance with 95%, 90%, . . . , 5% contours from outside to inside. The inset shows 
the distribution marginalized over the 10-11 kpc range in Rcc- 



Finally, the bottom panel shows the radial velocity (RV) versus Rcc distribution from 
combining the PDFs in radial velocity and distance for each star. The uncertainties in RV 
are about 3.2 km s" 1 and in Rcc are about 14.6% and smooth the total distribution. Even 
so, the velocity distribution shows two significant peaks around 10-11 kpc at RVs of about 
+27kms~ 1 and — 4kms~ 1 . 
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Discussion 



3.1. Dusty arms 



The jump in extinction around 9.5 kpc in the top panel of Fig. [2] indicates a broad and 
dense dust regi on in the path of the line-of-sight. Its location is just short of the Perseus 
arm at 9.9 kpc (IReid et al. 1 120091 ) implying that it is the dustlane of the spiral structure. 
Indeed, taking into account the Galactic rotation direction, the dustlane is in front of the 
stellar arm. Similar ly, the jump in ex tinction around 12.5 kpc is in front of the Outer arm 
located at 13.1 kpc feeid et al. Il2009h . 



3.2. Arm-interarm contrast 



In the anti-center direction the measured radial velocity RV is approximately equal 
to the Galactic space or cartesian velocity U, but with opposite sign. The bump in cr RV 
around 10 kpc indicates the mass center of the Perseus arm. Since the velocity dispersion in 
a equilibrium system is related to the mass density, the arm-interarm density contrast can be 
estimated. The peak value of a^y is 29.4kms _1 at 10 kpc, while it decreases to 27.8 km s -1 
at the interarm radius of around 11 kpc. Assuming an isothermal density distribution for 
the stars and a constant vertical-to-radial velocity dispersion ratio cr z /(TR, <Jr is proportional 
to X 1 / 2 where X is the surface density. The arm-inte rarm contrast is then K ~ 1.30. 
Alternatively, if a R is proportional to X 1 / 3 (IDehnen Ill999l ). we obt ain K ~ 1.5. The resulting 
contrast of K ~ 1.3 — 1.5 is in agreement with previou s studies (IDrimmel fc Spergel 1 12001 



Benjamin et al. 1 12005c iGrosbfll. Patsis fc Pompei 1 12004] ) . 



3.3. Bifurcation due to a resonance? 

The bifurcation in RV around 10-11 kpc could have different causes: (1) local dynamics 
of the spiral structure, (2) resonance of the spiral structure, or (3) resonance of the bar 
potential. 



Quillen et al. I (120111 ) argued that, within a density wave, the orientations of the orbits 
of the stars on the spiral arm are correlated with the mean radii of the orbits and result 
in an arc-like structure in the U — V velocity distribution. Since what we see is only the 



1 Here K = (S + $£)/(£ — <5£), in which <5£ is the incremental surface density in the arm, following 
Antoia et all (|201ll ) 



-7- 



marginalized distribution over V, an arc would be difficult to show. Even if the arc is 
elongated along V and thus show a sharp peak in U, it is hard to explain the ~30kms _1 
gap between the two peaks. Moreover, the bifurcation is at a larger radius than the location 
of the Perseus arm. Therefore, it seems that the local orbital distribution in the spiral arm 
cannot be the main reason of the bifurcation. 

If 10-11 kpc corresponds to the corotation ra dius of the spira l structure t hen its pattern 



speed is 20-22 km s 1 kpc 1 , in agreement with iGerhard I (120111 ). However, lAntoja et al. 



(120111 ) show that spiral structure alone results typically in too little substructure in the 
velocity distribution. 

The bifurcation can be the result of a resonance with the Galactic bar. The pattern speed 
of the Galactic bar has been reported in the range from 30 to 65 km s _1 kpc -1 depending on 
different data and method s (IGerhard 1120111). If we adopt 40 km s~ l kpc -1 and a flat rotation 
curve of v c = 220 km s -1 ([Rodriguez- Fernandez fc Combes 1 120081 ). the OLR radius is right 
around 10-11 kpc. It is then expected that stars group into inward and outward streams 
under the OLR radius, naturally explaining the bifurcation in RV. Meanwhile, the CR 
radius corresponding to this pattern speed is located around 5-6 kpc, which can explain the 
Hercules stream in the Solar neighborhood. As even the bar pattern speed itself is not yet 
well determined, further investigation in both observations and theory are needed to verify 
the resonant origin of the bifurcation. 



4. Summary 

In this letter, we use 698 red clump stars in the mid-plane of the Milky Way to derive 
the radial profile of extinction and radial velocity dispersion, as well as the velocity-distance 
distribution in the anti-center direction. We find: 

i) Dustlanes in front of the Perseus and Outer arms, consistent with external spiral 
galaxies; 

ii) A Perseus arm-interarm contrast between 1.3 and 1.5, in line with previous studies; 

iii) A strong bifurcation in the radial velocity distribution around 10-11 kpc, likely 
related to the outer Lindblad resonance of the Galactic bar, although additional observations 
are needed to constrain a possible spiral corotation resonance. 
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